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Abstract 

(— ( 1 Topological singularities occur in a broad range of physical systems, including collapsing 

stars and pinching fluid interfaces. They are important for being able to concentrate energy 
into a small region. Underwater air bubbles in particular appear in many practical applications, 
including new technologies to reduce skin drag on cargo ships. Previous theories show that 
just before an air bubble pinches off, the neck looks like a cylinder at its very smallest point. 
Unusually, however, the neck approaches this shape so gradually that the theoretical cylinder 
solution is not reached in practice; the singularity spends its entire lifetime in a transient phase. 
Therefore, in order to understand the evolution, we study the transient effects in detail. This 
CZ3 " paper details the simulation results of bubbles with initial conditions far from the cylindrical 

solution: squat, up-down asymmetric neck shapes, with imposed vertical flow. We hnd that 
the asymmetry is transient: the neck quickly shifts vertically to become up-down symmetric. 
Importantly, we find that the resulting symmetric singularity is a blend of the initial top and 
bottom sides, with a weighting factor that is tunable by adjusting the airflow through the neck. 
This effect should have implications for the later stages of evolution, including the generation 
of satellite bubbles and the formation of the Worthington jet. 

1 Introduction 

Non- linear continuous media can produce finite-time singularities, where a physical quantity nearly 
blows up in a fixed amount of time. In practice, the system never becomes fully singular; new forces 
become important at small scales and cut off the singularity. Examples of physical singularities 
include black hole formation, supernova explosion, crack propagation, and fluid pinch-off. Here, 
we study a singularity that occurs when an underwater air bubble pinches off from its air source. 
Despite the small energies i nvolved, the pinch-off produc es a small region of water moving at speeds 



Or 



produce 

of 10s of meters per second (jPlesset fc Prosperett? . 19771) . Many of the interesting features of bubble 



pinch-off apply to implosion- type singularities in general. 



Bubble pinch-off also has practical applications itself. For example. iMcCormick fc Bhattacharvva 



1973h discovered that injecting air bubbles under the hull of a ship reduces skin drag. Only re- 



cently has the idea been implemented on commercial ships, resulting in an overall 5-10% fuel savings 
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Figure 1: Bubble (dark) pinching off from a nozzle submerged in water (light). Bright spots are 
optical artefacts, (a) Initially, the bubble neck shape near the minimum has a generic quadratic 
profile. Zooming in (b-c), just before pinch-off the shape becomes two cones connected at the vertex 
by a short segment. A characteristic vertical length scale L z for the neck is the distance from the 
minimum to a height where the neck radius is 2i? m ; n (b) . Even at microscopic scales (c) the neck is 
relatively squat (L z /R m i n m 3.2, 9 + « 28°), still far from the theoretical end-state (L z /R m i n —> oo 
and 9 — > 0°). Neck vibrations eventually cut off the dynamics before the singularity is reached. 
(Images courtesy of N. C. Keim and S. R. Nagel). 



I Ferrante fe Elghobash i. 2004; Kawabuchi et «Ll . l201ll ). The drag reduction depends on the size dis 



tribution of the bubbles, which i n turn depends on the specifics of pinch-off (jLonguet-Higgins et al. 



199lHGordillo fe Fontelosl 12007!) . Al so, the bubbles can pinch apart in a turbulent flow, generating 



acoustic emissions (jFrizell fc Arndtl . 1987 ) and dispersing into a distribution of smaller bubbles 



which depends on the details of pinch-off ( Lasheras et al. . 2002). The mechanism s controllin; 



bubb le pinch-off are size-independent and also govern the breach of a submarine hull (llkeda et al. 
2009), and the violent whoosh behind an object plunged into water ijBergmann et all . 120061 ) . This 
paper focuses specifically on the effects of up-down asymmetry in bubble dynamics. Up-down 
asymmetry is present in nearly all bubbles due to the hydrostatic pressure gradient. In general, the 
bottom portion of the neck is more slender than the top due to the hydrostatic pressure. This effect 
can b e exacerbated if the bubbles are generated in a dynamic stream of water (jOguz fc Prosperettil . 
1993). In fact, this research has potential implicat ions for the jet of wat e r that s hoots upward after 



such an impact, called the Worthington jet after IWorthington fc Cold (|l897l ). (jGekle fc Gordilld . 
20ioh , S3 discusses how the initial conditions of the cavity determine the bubble shape just before 



the impact that spawns the Worthington jet. 

Few people realise that they are creating a singularity when they blow a bubble into a glass 
of water using a drinking straw. In experimental set-ups, the effect can be controlled by precisely 
pumping air into a pressurised nozzle (see figure [TJl) . As the air rises, surface tension competes with 
buoyancy to create a quadratic neck shape. Then, the neck becomes unstable and water rushes in 
from all sides to close it off. This sets up the regime we study in this paper, where the inertial flows 
diverge in speed where they are closing the shrinking neck, dominating all other forces and creating 
a singularity. If this regime were allowed to co ntinue unaltered, the neck wo uld approach a long- 
and-slender cylinder, with purely radial flows (jLonguet-Higgins et all . 119911 ). Bef ore this occu rs, 



however, fast air-flows develop in the neck, potentially altering the collaps e scali ng (Gordillo et al.. . 
2003lBurton fc TaborekL l2008h . Eventually, the air even g oes supersonic jGekle aZ.l.l2010h . Just 



before the airflow effect, some authors including iGordillo et all (|2005h and iKeiirJ j 2 0111 ) measure 
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a small residual asymmetry that influences the air flow. This is consistent with our results that 
the initi al asymmetry deca ys as a power law, but is still measurably appreciable after at this scale. 
Finally, iKeim et all I^POrt )" found that if the bubble is not quite axisymmetric to begin with, the 
deviation actually becomes more pronounced as the neck approaches pinch-off. The outcome is a 
side-to-side contact before the average radius reach es 0, cutting off the si ngularity. The se vibrations 
were confirmed and st udied in detail in l ater theory dSchmidt et al. , 2009h . experiments ( Keim . 2011 ; 
Enriquez et al. , 2011 ) , and simulations ( Turitsvn et al. , 20091 : Lai , 2012f) . 

This singularity is unusual in that it reaches the cut-off point before the neck shape has a chance 
to reach its predicted end-state shape, which is a long cylinder (connecting slender cones on the 
top and bottom, see §2). Most singularities with a unique e nd-state shape ( universal singularities) 
rapidly approach that shape with power- law dependence ( Shi et all . Il994h . In those cases, the 
transient effects leading up to the universal shape are short-lived and ignorable. In contrast, the 
air bubble neck shape evolves very gradually as a function of the nec k radius i? m in, only going as 
log (i? m in) (|Gordillo et all 120051 : lEggers et all l2007t iGekle et all l2009h . Indeed, the transition is so 
slow that on observable length scales the neck never appears cylindrical, but rather as two relatively 
squat cones, smoothly connected (see figure [T|>c). Although technically universal, this singularity 
spends its lifetime in the transient phase. Therefore, this paper explores the various transients that 
can occur in evolution. 

We begin by exploring the effect of highly distorted initial neck shapes, including necks with large 
curvature and up-down asymmetry. We found that the asymmetry is truly transient in every sense: 
it vanishes after only several factors of 10 decrease in the minimum radius. The neck undergoes 
this transition not by flexing the cones, but by shifting the vertical position of the pinch-off. By 
shifting instead of readjusting in place, the (new) collapse region rapidly becomes nearly up-down 
symmetric. Then, the neck goes through a phase of being symmetric but also squat, with significant 
(symmetric) vertical flows ( Herbst fc Zhangl . 1 2 llh . Depending on where the neck shifts vertically, 
the symmetric squat neck inherits its dynamics from either its top or bottom side, or a blend of both. 
Interestingly, we found that by tuning the initial vertical flow near the neck, we could adjust the 
degree of this blending. This effect ma y explain some of the late-stage effects seen in experiments 
that tune the density of the inner gas ( Burton fe Taborek . 2008t Keim . 2011 ). When the bubble 
is in the formation stage, gas rises through the neck. Denser gasses impart more momentum to 
the neck area, giving the neck an initial vertical velocity component. We pred ict that this should 
cause the neck to inherit its aspect ratio from the top side, which is more squat. iBurton fc Taborek 
(2008) do in fact see more squ at end-state neck shapes for higher gas densities (up to their observed 
critical density). IKeim ( 2011 ) observes that the very early stage air flow during formation sets up 
a water velocity field, which then evolves independent of the air during the inertial stage studied 
in this paper. He also confirms our observations that the water flow affects the up-down symmetry 
and neck position later in the evolution. Neither experiment, however, attempted to manipulate 
the initial neck velocity field purposely, except by altering the gas density. Given our simulation 
results, it should be worthwhile to explore more deliberate ways of controlling the initial airflow 
through the neck, for example by blowing the bubble with a top and bottom nozzle simultaneously. 

^2] summarises the previous research on the inertial stage of the bubble singularity, and the need 
for looking at up-down asymmetry separately. $3] gives an overview of the computational methods 
used to perform the studies. The results are divided into two sections for ease of reading: §31 gives a 
full description of a typical collapse, ^describes how the initial flow field, specifically the vertical 
velocity gradient, can be used to tune the details of the final collapse. i}5] concludes with a discussion 
of how the results could be seen in experiment, and problems that remain unsolved. 
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2 Background 



The first stage, described in SJTJ se ts up the initial conditions for the implosion, and involves the 
interplay between surface tension and buoyancy. After this point, surface tension becomes dwarfed 
by inertial forces. This is the point where we begin our simulations, and we assume inviscid flow 
with only inertial and pressure forces. 

We use the common assumption that the water's velocity field u is incompressible (V • u = 0), 
irrotational (V x u = 0), and decays to zero far away. 

We add another assumption that the air in the bubble is dynamically passive, with a uniform 
pressure P a ir(i) whose value ensures constant bubble volume in time. Experiments show that after 
the neck begins to collapse, airflow is u nimportant until the neck radius reaches roughly 100 pm 



(jBurton fc Taborekl . 12008c iKeiml . 1201 lh . Around that point, airflow through the neck becomes 
relevant and our assumption breaks down. Soon after, neck vibrations dominate. Here, we will 
only study what happens up to this point, assuming an axisymmetric bubble with uniform air 
pressure. 

Being curl-free, the velocity can be described by a scalar potential, defined as V</> = u. Since 
the exterior flow is incompressible, V- (V(/>) = implies 



V z d> = 0. 



(1) 



That is, Laplace's equation governs the potential in the bulk. This allows u to be solved everywhere 
in the exterior given <f> on the interface. Equating the relevant stresses on the interface gives a non- 
linear differential equation for <f) on the surface of the bubble, first-order in time and space: 



0. 



(2) 



where p is the water density. The previous two equations determine the evolution of the velocity 
field. The velocity field in turn determines how the interface evolves. Surface elements are advected 
by the local velocity on the surface, S : 



dx 
ch 



= v« 



is • 



(3) 



Previously, the sol utions to equat i on [21 were only studied carefully in the up- down symmetric, 
long-and-slender limit (jGordillo et al , 2005 ; Eggers et ai , 2007 ; Gekle et aUl2009h . These solutions 
produce elegant, analytic results showing the neck approaching a cylinder logarithmically slowly as 
-Rmin — > . There are differences between the theories, but each assumes a (time-dependent) purely 
radial flow. This is equivalent to saying that the potential distribution is proportional to the zeroth 
order cylindrical harmonic: 



ficyi (*) = -q (t) log [r/r-oo (t)] , 

u cyl (t) = -q(t)f/r, (4) 

where r M (t) is a length scale related to the system size, and q (t) is related to the rate of inflow, and 
is determined by equation [2] For a perfectly cylindrical neck with constant P a j r and , the reader 

may verify that q(t) = [(P a i r i?^ nin + const) / log (r - 00 /i? m i n )] 2 satisfies the governing equations. 
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Figure 2: The near-cylindrical theory of Gekle et a l. (2009) (dotted lines) fails to properly reproduce 
the simulation (solid lines, see §3]) for a neck shape that is initially up-down asymmetric and squat 
(6 + = 40° and 6~ = 5°), with initial radial flow. The simulation reveals that the asymmetry is a 
short-lived transient, whereas the theory predicts that the top and bottom sides evolve slowly and 
separately to 0° cone angle. 



In most cases, the pressure term is negligible, leaving q(t) ~ [log (roo/Rmin)] 2 , which is nearly 
constant. This defines a flow, u = — [log (r^/R ^r,)] -2 r" 1 ^, that evolves th e neck gradually 



111112 
Longuet-Higgins et al\ (|l99ll ) found that any 



parabola-of-rotation would generically evolve into a hyperbola-of-rotation. Indeed, at late stages 
the bubble neck shape does resemble a hyperbola-of-rotation by the fact that the neck looks like 
two smoothly connected cones (figure [TJ;), which we classify by top and bottom opening angles, + 
and 6~ respectively. For consistency, we actually measure an effective cone angle 8 e s defined by 

2-Rmin = Rmin "J 1 + (tan 2 6 e e) (L z /Rmin) 2 ■ Taking into account the logarithmic dependence of q on 
i? m i n , the cone opening angles decrease very slowly to 0, approaching a cylinder. Even at the small 
scale of figure [TJ;, the cone angle is approximately 28°. Even under ideal conditions the neck does 
not become cylindrical before neck vibrations tear off the bubble. 

There are two problems with analytic solutions of the form in equation 2] First, they only apply 
when the neck shape is long-and-slender and up-down symmetric. Neither assumption holds for real 
bubble necks (see figure [1} . The theory predicts qualitatively incorrect behaviour for initially up- 
down asymmetric necks (figure^. In the theory, the neck shape evolves as a stack of non- interacting 
cross sections, each one governed by equation [4] (or similar). In reality, vertical interaction plays 
a large role in the shape evolution. Without vertical flow, the asymmetry decays logarithmically 
slowly. With vertical flow, we find that the asymmetry vanishes as a power law (measured to go as 

n-l \ 
-^min/ - 

In fact, for up-down asymmetric shapes, the most important effect comes from the global 
interaction between the top and bottom sides. If we artificially forbid flow through the mid-plane, 
z = 0, the up-down asymmetry no longer vanishes quickly (figure [3] dashed lines). §31 describes the 
mechanism in detail. In short, the singularity shifts vertically in the direction of the smaller cone 
angle, and in doing so equalises the cone angles. §31 shows that this effect can be manipulated. By 
introducing a vertical flow to the initial conditions, the vertical position of the neck can be tuned. 
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Figure 3: Only when the top and bottom sides fully interact does our simulation show the up-down 
asymmetry to be short-lived (solid lines). For comparison, we also evolved two symmetric shapes 
with + = 9~ = 40° and 6 + = 9~ = 5°, respectively (dashed lines). Taken together, these two 
simulations show how the asymmetric shape would evolve if the vertical position of the neck were 
fixed (and no flow were allowed through the z = plane). In that case, the asymmetry decreases 
very gradually. This comparison shows that the asymmetry vanishes quickly because the top and 
bottom interact. Specifically, the vertical position of the neck shifts. 



In doing so, the resulting symmetric neck inherits a different blend of characteristics from the top 
and bottom sides. This effect should be realisable in experiments by adjusting the air flow rate 
through the neck during the initial set-up phase. 



3 Methods 

To Solve equations [TJ(3j we use the standard procedure of reformulating the equations as a boundary 
value problem, where the bulk differential equation becomes a surface integral equation. This re- 
quires two steps. First, the ordinary derivative in equation[5]is replaced by the co-moving derivative 
Tn ~ Th JrU ' ana ~ w ^ n a PPly to the fluid parcels that constitute the surface. Second, equation 
[J is reformulated using Green's integral: 

,/ s I ( (xi -x o )0(xi) VjHxi) 

27T</>(X )= O) ■ -3 + 1 

J \ |xi-x | d |xi-x | 
xies 

where xo is a fixed point on S, xi parameterizes the surface, and n is the unit surface normal at 
xi pointing into the bubble. 

With both governing equations defined on the interface, the formulation reduces to 2 dimensions. 
We further reduce to I dimension by assuming axisymmetry. At each time step, we use cubic 
splines to interpolate the N(t) discrete surface nodes Xj. Then, for each i 6 {I . . .TV} we allow 
Xj to take the place of xo in equation ([5]). We use Gaussian quadrature to perform the integral 
over segments connecting adjacent spline midpoints j € {I . . . N} , with the first and last segments 
being only half splines. Thus, the N instances of equation [S] can be schematically written as 
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Figure 4: (a) The initial surface is comprised of an upper half with opening angle 9 + — 40°, and a 
lower half with opening angle 6" = 5°. The initial flow is radial, (b) Zooming in 50 x , four snapshots 
show the neck becoming rapidly more symmetric via a descent of the neck position (indicated by 
black curves). In this egregiously asymmetric example, a second neck actually forms causing a jump 
in the vertical neck position (seen as a tiny kink in the black curve in b). (c) Zooming in another 
50 x, the simulation evolves the bubble to length scales smaller than physically possible. Here, the 
shape is nearly cylindrical and the neck shift is insignificant. Time steps between curves are: (b) 
0.19i?o/wor: an d (c) 4.4 • 10~ 4 R /u 0r . Each innermost profile becomes the outermost in the next 
image. Every tenth node is shown. 



fa = Aij(f)j + Bij (it_i_)j , where A and B are N-by-N matrices and u± = • n. Therefore, given 
<f) on the surface we have u± — B^ 1 (fa — Aij(f)j) and u\\ = where ds is an infinitesimal arc 
length. Using At ~ Rm\n/Rm\n , we advect the N nodes using u and evolve cf> using equation [2j 
completing the cycle ( Pozrikidis . 1997 : Oguz fc Prosperettil . 1993f ). The pressure P a j r can be solved 



for explicitly after evolving the shape and before updating fa (formula and derivation available 
upon request). 

In order to accurately resolve the pinch-off dynamics, we found it important to use a node 
distribution scheme that maintains a gradual variation in the spacing between node points and 
that continually adds nodes in the neighbourhood of the minimum (see figure |4j. At each moment, 
if the vertical distance Az between nodes exceeds a maximum spacing Az max (i, z) in the region 
where the neck radius is less than 4i? m i n , we add a new node at the midpoint of that spline. After 
experimenting with several functions for Az max , we found Az max = max [0.15 min (|z| , L z ) , 0.08L Z ] 
allows the simulation to accurately track the dynamics 0. 

We begin the simulation by specifying an interface shape and velocity field. In general there 
are an infinite number of each. Here, we examine generic smooth distributions. To expedite the 
computation, we do not prescribe a purely quadratic shape profile at t = 0, but instead use the 
result from previous studies that a slender quadratic neck evolves into a hyperbola, and prescribe: 



R(z) = J(Rl + zHan 2 9) - (z 4 /h 4 ) {Rl + h 2 t&v 2 9), 



1 Fov the up-down asymmetric case, we define the current height of the neck minimum as zo, R (^zq ± L^j = 
2/? m ; n , and A2 max = max |^0.15 min (\z — Zq \ , , 0.08 min J, L^~^ J for z > zq (reverse superscripts for z < zq). 
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where 29 is the opening angle of the cone, Rq = R m i n (t — 0) , and h = 140i?o is height from the 
neck minimum to the end of the bubble. 

To initialise the velocity field, most previous research uses a pure radial collapse, because it is 
the counterpart of the cylindrical end-state solution. We use this as a starting point. Equivalently, 
we prescribe a line of point sinks along the z-axis: 



$U. (r,z) 



r 2 + (z-(Y 



dC 



— L 



In 



Z + L + Jr 2 + + L) 



z- L + yr 2 + (z- L) 

We use L = 100, much larger than the neck region, to avoid end effects. For we also add a 
vertical velocity component to the initial conditions of the neck region. One trivial possibility is a 
uniform vertical component. That would be the same, however, as just boosting the neck region to 
a new reference frame z = UQ Z t, and the pinch-off would be identical, only boosted. Instead, we use 
the following form, which consists of a line of sinks (sources) along the z-axis that vary in strength 
proportional to height £ : 



r — — 

Mr,*) = / C(r 2 + (z~0 2 ) "dC 



aL — d + + z log [L — z + d-] — z log [— L — z + d + ] , 



(G) 



where d± = y r 2 + (L ± z) , and the limit r — > is taken for points on the z-axis. This field 
acts as a stretched out dipole. The total initial potential field consists of a weighted sum of the line 
sink and the dipole field: 



®tot = 



Rodo 
2L 



UQr^l.s. 



2L , 
■— +log 
do 



u 0z <S> z 



(7) 



where do = \/ Rq + L 2 . The terms in parentheses are normalisation factors, so that the initial 
velocity at the neck minimum (z = 0) has radial and vertical components u 0r and Uq Z7 respectively. 
This initial velocity field has the added benefit of resembling the effect of vertical airflow through 
the neck on the surrounding water. 



4 Results 

Our simulations show that the evolution of the neck towards singularity depends critically on both 
the initial neck shape and the initial velocity field. We show the two effects separately. 

21 details how an up-down asymmetric shape evolves with a simple initial velocity field: pure 
radial flow. The velocity field immediately develops a pressure ring just outside the neck region, 
but below the neck minimum itself (in general, towards the side with smaller cone angle). As the 
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Figure 5: The characteristic pressure peak that drives the singularity develops below neck minimum 
of the initial shape (original neck shape: dark dashed curves; original neck minimum: thin horizontal 
dotted line), (left). The pressure peak lies just outside the new neck minimum, with steep gradients 
between the peak and the surface, along which the pressure is fixed to be P a ; r (t) , (right). (Note: 
this point of evolution corresponds to the centre- most profile in figure HJd) . 

pressure peak moves inward, it bulges the shape inward, shifting the vertical location of the neck 
minimum and making the neck symmetric. The symmetric neck then progresses as a symmetric 
collapse, with the upper cone angle adjusted to match the bottom. 

33] adds the effect of initial vertical flow around the neck. A small upward flow is enough to 
oppose the neck from shifting downward. After becoming symmetric, the neck is a blend of the top 
and bottom cone angles. This gives a possible experimental method for adjusting the cone angle of 
the neck at pinch-off. 

Shape Asymmetry 

We studied many different combinations of initial top and bottom cone angles (8 + and respect- 
ively), with each angle ranging between 5° and 60°, along with different initial velocity distributions. 
The results are qualitatively the same for any combination of top and bottom cone angles, but the 
asymmetric effects become more pronounced for larger asymmetry. In this section, we use a neutral 
radial inflow (equation [7] with uo z = 0), and an initial 40° top cone and 5° bottom cone. 

The initial shape consists of the neck, a hyperbola-of-rotation with the given cone angles, roun- 
ded off by large end-caps (figure[2Ji). Immediately after the simulation starts, a ring of high pressure 
develops in the water surrounding the neck, but slightly below the initial neck minimum. This causes 
the shape to bulge in at a height below the original neck minimum (figure HJd), and the position of 
the neck minimum to descend. In our extreme example, the descent occurs so fast that a new neck 
actually forms, which appears as a jump in the neck's vertical position (figure HJd). Meanwhile, 
the pressure peak continues to grow in strength and become more symmetric (figure [5]) • The neck 
minimum shifts downward at roughly a constant rate (lab frame) until pinch-off (figure HJd). This 
constant vertical rate is soon overcome by the large radial implosion flow, and becomes insignificant 
at later times (figure [4]:). By this time, the neighbourhood of the neck is completely symmetric, 
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Figure 6: (Left) The pressure peak from figure [5] creates a symmetric stretching flow centred below 
the original neck minimum, with positive vertical velocity component u • z above the peak and 
negative below. This flow leads the shape into a symmetric, vertically-offset collapse. (Right) 
Initially, the vertical flow is prescribed with no vertical flow near the neck minimum (innermost 1). 
Immediately, the pressure peak (figure [5]) develops, mainly generating positive vertical velocities 
on the top side (2). The pressure peak shifts downward, forming a new, symmetric neck, below 
the original minimum (3). Out of this neck, an up-down symmetric straining flow develops (4-7). 
Curve 5 (thick) corresponds to the left figure. When the straining flow begins to develop, it first 
is only a local maximum on the top surface (*), but continues to grow (6-7), eclipsing the original 
peak on the top side, which stagnates and becomes irrelevant. The straining flow is slightly offset 
in u z , signifying that the neck continues shifting at a constant rate in the lab frame. 

and the system undergoes the previously-studied symmetric collapse. In this process, the top side 
of the shape performs most of the accommodation to match the bottom side, and the large initial 
upper angle has little effect after the transient (figure \S§ . 

The process can be seen in more detail by studying the evolution of the velocity field, specifically 
the vertical velocity. Initially, no vertical flow is present (figure [6] curve 1), but immediately an 
asymmetric straining flow develops (figure |6] curve 2). Soon afterwords, the origin of the straining 
flow shifts downward (figure [6] curve 3). Driven by the pressure ring, the straining flow becomes 
symmetric (figure [6] curve 4), which in turn makes the neck shape symmetric. Eventually, the 
vertical velocity forms a local maximum in u z on the upper side (figure [S] curve 5, indicated by 
'*'). From then on, the straining flow grows asymptotically strong, leaving behind only vestiges of 
the original top side, which stagnate and become insignificant (figure |6] curves 6-7). The transition 
from curve 2 to curve 4 shows that the top part of the straining flow adjusts to match the bottom 
half through a downward shift of the straining flow centre. 

The shape and velocity field are interdependent, and it is informative to show how they co-evolve. 
Figure [7] condenses this information into a phase space. The vertical axis gives the ratio of the 
maximum surface vertical velocity component to the maximum surface radial velocity component 
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symmetric shapes: 




Figure 7: The two sides of the asymmetric evolution (solid curves) evolve drastically differently 
from their symmetric counterparts when viewed in phase space (velocity ratio (max \ u z \ /max |w r |) > 
versus the neck aspect ratio L z /R m - m ). Due to the strong asymmetry, the two sides see each other 
and their phase trajectories attract. The 40° side is stretched as the neck rolls down, while the 
5° side is squashed. Meanwhile, the 40° side is prevented from reaching its full potential vertical 
velocity. This is because the singular straining flow is generated on the 5° side (see figure [5]). The 
newly formed local maximum vertical velocity (figure [6J '*') is represented here also by a '*'. This 
peak eventually prevails as the global maximum. The two phase curves finally meet, as the neck 
region becomes symmetric, then slowly becomes more slender going towards pinch-off. 



(two distinct points). The horizontal axis gives the shape aspect ratio L z /R m { n , i.e., when the 
trajectory moves to the right, it corresponds to the neck shape becoming more long-and-slender. 
The top and bottom sides each have separate trajectories, and for comparison, we show each side's 
symmetric counterpart (which shows what would happen if no flow were allowed through the plane 
z = 0). Each trajectory starts on the horizontal axis, corresponding to the prescribed zero initial 
vertical velocity, at a position depending on the initial shape. For example, the 40° initialised side 
begins with a corresponding aspect ratio of L z /Rp = y/3 cot 40° = 2.0 6. After an initial transient, 
all the curves coalesce onto a universal curve (jHerbst &: Zhana 120111 ) and slowly progress to the 
universal form: long-and-slender (large shape aspect ratio) and radial flow (zero vertical velocity). 
The progression along the universal curve has been studied. We focus here on the transient. 

Each side of the asymmetric shape differs significantly from its symmetric counterpart, indicating 
that the vertical flow through the mid-plane (interaction between the top and bottom) z = is 
important. The flow through the mid-plane allows the neck to shift until the straining flow is 
symmetric. The result is that the trajectories of the top and bottom sides "attract" one another in 
phase space. In this case, the direction of the neck shift is downward. The downward flow of the 
neck shift opposes the straining flow on the 40°-initial top side, and augments the straining flow on 
the 5°-initial bottom side. The result is that the peak surface vertical velocity ratio undershoots 
its symmetric counterpart on the top side (figure [7] top solid vs. dashed), whereas it overshoots its 
symmetric counterpart on the bottom side (figure [7] bottom solid vs. dashed). Similarly, the squat 
top side shape becomes more slender and the slender bottom side becomes more squat. In this way, 
the two trajectories are brought together in a short time. Some important features seen in previous 
plots can also be seen in figure[7J The local maximum vertical velocity ('*' in figure^ is shown as 
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10" 3 10" 2 10" 1 10° 

Figure 8: Adding some vertical flow (uq z = 0.4uor) to the initial conditions causes a significant 
change later in the evolution. Instead of following the trajectory of a 5°-initial shape (see figure [3]), 
the evolution more closely follows the 40°-initial shape. 

a separate trajectory in figure [7] (also indicated by '*'). Also, figure [7J shows that the bottom side 
does not have to detour much in order to meet up with the top side trajectory. This effect was also 
seen in figure [3] In short, since the neck shifts down in this case, the singularity mostly inherits the 
characteristics of the bottom side. It is the top side that adjusts to mirror the bottom side, and 
the singularity ignores the fact that the top was initially squat. This effect can be reversed, as seen 
in the next section. 

Tunable cone angle 

To understand the context for this section, refer back to figure[3J There, after a short transient, the 
top and bottom sides (solid lines) converge to a single trajectory that closely follows the trajectory 
of a 5°-symmetric initial shape (lower dashed line). Here, by adding a vertical velocity to the 
initial conditions (equation [Jj with non-zero uq z ), the picture changes. A prescribed initial vertical 
velocity 40% in strength to the initial radial velocity, uq z — 0Auo r , causes the shape to evolve into 
a symmetric shape that more closely follows the top-side 40°-symmetric initial trajectory (figure [51 
compare to figure [3]) . The added vertical velocity opposes the tendency for the neck to shift down 
in the early evolution. As a result, when the neck becomes symmetric, it inherits the characteristics 
of the top side, rather than the bottom side. This modest change in initial condition causes a 
qualitatively different outcome. 

The effect is not bimodal, but continuous (see figure^. For UQ Z /uo r between the range of about 
0.0 to 0.5, the resulting trajectory smoothly varies between the bottom-side dynamics and the top- 
side dynamics. For initial downward flow {uq z < 0.0), the system is insensitive to the exact value; 
any downward velocity just results in the singularity inheriting the bottom-side dynamics (there is 
no overshoot effect). Interestingly, the dynamics never fully approach the top-side trajectory, but 
seem to asymptote slightly below it. We do not have an explanation for this disparity. 

With an added vertical velocity UQ Z /uo r — 0.4, it is the top side trajectory that dictates the 
evolution, with the bottom side going out of its way to accommodate (figure [TU]). The point on the 
phase diagram where the two trajectories merge (becoming symmetric) is shifted to the left and up, 
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Figure 9: For neutral and downward initial neck velocities, the resulting shape after the transient 
(fixed here as -R m ; n = 0.001i?o) is the same as if the initial shape were 5°-5° symmetric (lower 
line). That is, the upper, 40° side has no effect on the late-stage dynamics. For upward initial neck 
velocities, i.e. uq z > 0, a new effect occurs. The upward flow disrupts the downward neck shift, and 
the post-transient neck shape is a blend of the initial top and bottom sides. Evidently, a modest 
initial vertical flow around the neck can tune the properties of the end-state singularity. 




Figure 10: Adding a small vertical flow to the initial conditions has a big consequence on the 
transient. Here, the initial shape is the same 6+ = 40°, = 5°, but now there is an added 
vertical flow given by equation [7] with itoz/wor = 0.4. The main effect is that the point where the 
two solid lines meet (the point where the shape becomes symmetric) is significantly further left 
(corresponding to a larger effective cone angle). The upward initial velocity opposes the downward 
neck shift, and the resulting singularity retains more of the squatness of the upper side. This effect 
is tunable by adjusting uo z /uor- 
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corresponding to a more squat aspect ratio with larger straining flow. Without the downward neck 
shift, the top side evolves almost as if it does not see the bottom (figure [TOl similarity between top 
solid and top dashed). The initial vertical flow keeps the neck position elevated, and the singularity 
develops within the context of the top side. The resulting singularity inherits its characteristics 
from the top side, nearly ignoring the bottom side completely. A small initial vertical velocity near 
the neck can change whether the singularity inherits from the top side, the bottom side, or a blend 
of each. 



5 Conclusion 

In general, bubble pinch-off begins with some up-down asymmetry. We found that the neck adjusts 
the vertical position of the pinch-off point to make the neck shape and velocity field symmetric. If 
the initial flow field around the neck is purely radial, then the neck shifts (down) towards the more 
slender cone. The new neck forms with a top that emulates the bottom. The resulting singularity 
depends only on the initial 0~ . The situation is reversed if a strong upward velocity is introduced to 
the neck during the creation stage. Then, the vertical flow prevents the neck from descending, and 
the new neck depends only on 9 + . For intermediate upward velocities, the shape inherits a blend 
of the two angles. 

If the singularity were allowed to proceed all the way to i? m i n = 0, this effect would not matter 



since all symmetric necks, even squat ones, progress to a cylinder eventually (jHerbst &i Zhang 



2011). In reality, however, the asymptotic form progresses so slowly that other dynamics pre-empt 
the singularity before the neck can become slender. Therefore, it is the post-transient shape that 
sets up the post-pinch-off effects. For example, we expect that tuning the post-transient cone angle 
will play a significant role in the later effects, including neck air flow, azimuthal vibrations, and the 
Worthington jet, but more research is needed. 

This work was supported by NSF No. CBET-0967282 and the Keck initiative for ultra-fast 
imaging (University of Chicago). The author thanks Justin Burton, Michelle Driscoll, Nathan Keim, 
Lipeng Lai, Sidney Nagel, Thomas Witten, and Wendy Zhang for discussions and suggestions. 
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